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Abstract. - The properties of perfect screw dislocations have been investigated for several 
zinc-blende materials such as diamond, Si, /3-SiC, Ge and GaAs, by performing first principles 
calculations. For almost all elements, a core configuration belonging to shuffle set planes is fa- 
vored, in agreement with low temperature experiments. Only for diamond, a glide configuration 
has the lowest defect energy, thanks to an sp 2 hybridization in the core. 



Several technologically interesting covalent materials are semiconductors or insulators with 
a zinc-blende crystalline structure. The most famous is silicon, extensively used in electronic 
devices. It is usually considered as a model for other materials with the same structure, and 
most of its properties are well known, thanks to an impressive number of dedicated studies. 
For instance, the plasticity of silicon has been largely investigated both experimentally and 
theoretically. At high temperature, silicon is ductile, and dislocations are dissociated into 
Shockley partials; since the dislocations are frequently aligned along the (110) dense directions, 
they mostly appear as 30° and 90° partials and slip in the so-called 'glide' set of {111} planes 
(Fig. [T|) [1]. Calculations have allowed to determine the core structure of partials [2-4], 
although the most stable configuration of the 90° is still not known with certainty [5-8]. At 
low temperature, in the brittle regime, surface scratch tests and confining pressure experiments 
have recently revealed the presence of undissociated dislocations, with screw, 60°, 30° and 41° 
orientations [9, 10], so that very little investigation of these perfect dislocations recently have 
been reported. So far, theoretical studies have shown that the non-dissociated screw belongs to 
the 'shuffle' set, though there has been a controversy about the most stable structure [11,12]. 

Other covalent materials with zinc-blende structure are expected to behave like silicon, 
at least at high temperature. Partial dislocations have been theoretically characterized for 
/3-SiC (cubic 3C-SiC) [13-15], diamond [16-18], Ge [7,8,19], GaAs [20,21], and others [22,23]. 
However, it is not clear whether the results obtained at low temperature for silicon remained 
valid for those materials, and little is known on that subject. Experiments have shown that 
low temperature deformation of GaAs yields undissociated screw dislocations [24], like in 
silicon, and recent calculations by Blumenau et al focussed on perfect screw in diamond [15]. 
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Fig. 1 - Ball-and-stick representation of the zinc-blende structure ((101) plane). The A, B and C 
circles indicate three possible positions of the dislocation line. Broken (dotted) lines show the 'shuffle' 
('glide') {111} planes. 

Recently, we have investigated the properties of undissociated screw dislocations in several 
zinc-blende materials [25]. Our simulations, performed with empirical potentials, indicated 
stable configurations in the shuffle set for Si and Ge, and in the glide set for /3-SiC and 
diamond. One has to be cautious with these results though, since empirical potentials may 
fail to describe correctly the energetics of the highly distorted core geometry [12]. This is 
especially true in the case of SiC or GaAs, for which charge transfers occur. 

In this paper, we report first principles calculations of undissociated screw dislocations 
in several zinc-blende materials, namely Ge, Si, /3-SiC, diamond and GaAs. The theoretical 
approach is very similar to the one proposed in [12]: a quadrupolar arrangement of screw 
dislocations is built using anisotropic elasticity theory, from which a 12 x 12 x 1 cell includ- 
ing two dislocations with opposite Burgers vectors is cut. Periodic boundary conditions are 
applied to recover the quadrupolar arrangement [3]. After relaxation, quantities attached to 
a single dislocation are obtained by removing contributions from dislocation interactions, de- 
termined from anisotropic elasticity theory [12]. Density functional theory calculations have 
been performed with the ABINIT code [26]. Norm-conserving pseudopotentials [27], local 
density approximation, and two special k-points along the dislocation line have been used. 
We employed plane-waves energy cut-off of 10 Ry for Si, Ge and GaAs, and 40 Ry for SiC 
and diamond. These parameters allow a satisfactory determination of elastic constants in all 
materials (Tab. IJ). 

Three possible dislocation line positions have been considered in this work (Fig. [TJ| . A is in 
the center of the projected hexagon and is at the intersection of two shuffle planes. Position 
B is at the center of a long hexagon bond and belongs to a shuffle and a glide planes at once. 
On the basis of geometrical arguments, it has been suggested that this mixed shuffle-glide 
configuration is a saddle point, and therefore, is not stable [28]. This has been confirmed for 
Si by first principles calculations [12]. Therefore, in this work, the B configuration is obtained 
by constraining the two atoms closest to the dislocation line to remain at the same distance 
along [101] . Finally, C is at the center of a short hexagon bond, at the intersection of two glide 
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Table I - Elastic constants, in Mbar, calculated for several materials and compared to experiments 
(in parentheses). 
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C44 


5.92 (5.78) 
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planes. The table [TT] shows different energetic parameters, deduced from the first principles 
calculations, such as the energy differences between configurations, the c ore r adii according 
to [1], and the core energies, obtained for a core radius equal to the Burgers vector. Data 
for silicon, already presented in [12], are included for comparison with other materials, and 
will not be discussed here at length. The most important point is the stability of the shuffle 
configuration A, both energy of B and C cores being substantially higher. 

In the case of diamond, it appears that the C core is clearly the most stable configuration, 
and B the less favored. The large energy difference with the A core may be explained by the 
narrowness of the C geometry, as well as the presence of two atoms in the core with a sp 2 
hybridization, energetically favored for diamond. The bond length between these two atoms 
is 1.31 A, less than the first neighbor distance in graphite (1.42 A). Our results indicate that 
a non dissociated screw should be located in glide planes. However, it has been shown that 
dissociation is favored in this set of planes [15], and non dissociated dislocations are then 
unlikely to happen in diamond. The determined energy ordering is in agreement with a recent 
study [15], but our computed energy differences are much lower, with AEac — 1-29 eV/Bv 
and AEbc — 2.51 eV/Bv to compare with 3.3 eV/Bv and 5.4 eV/Bv respectively. The lower 
precision of the self-consistent tight-binding method used in [15] may explain this difference. 

Silicon carbide is interesting because it can be considered as an intermediate case between 
silicon and diamond. Deformation at high stress/low temperature of hexagonal 4H-SiC sug- 
gested the coexistence of undissociated and partial dislocations [29]. However, apart from 
our previous semi-empirical approach, undissociated dislocations have not been investigated. 
Our calculations reveal that the A configuration is most stable like in silicon. The C con- 

Table II - Energetic data calculated for the A, B and C configurations: total energy in eV /Burgers 
vector (relative to the most stable configuration, in box), core radius, and core energy (calculated with 
a core radius equal to the Burgers vector). Note that the configuration C is not stable for Ge and 
GaAs and relaxes to A. 
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Fig. 2 - Variation of screw configuration energies, for Ge, Si, SiC and diamond, calculated with the 
Tersoff potential and with a first principles method. The A core energy is the reference for each 
material. 



figuration is second, with an energy difference much lower than for diamond. These results 
are in contrast with previous Tersoff potential calculations, giving C as the most stable core 
structure [25] . Screw configurations for silicon carbide include the same number of core C and 
Si atoms, due to the zinc-blende structure. Therefore, in the case of a C configuration, sp 2 
hybridized core atoms are one carbon and one silicon, the former being energetically favored 
unlike the latter. We have tried to build a C configuration with two sp 2 carbon atoms, by 
exchanging the center Si atom with one of its carbon first neighbor. However, the resulting 
structure is not favored, with an excess energy of 2.58 eV/Bv compared to the original C 
configuration. Non-stoichiometric core configurations have not been investigated, such a task 
being beyond the scope of the present study. 

Germanium and gallium arsenide show very similar behaviors regarding screw characteris- 
tics. In fact, only the A configuration is stable, and the energy differences A£ab are very close 
for both materials. Also, the A core was obtained after relaxation of an initial C configuration 
in both cases. Non dissociated screw dislocations have been observed in GaAs deformed at low 
temperature [24] . Our calculations unambiguously indicate that such dislocations are lying in 
shuffle set planes. 

In the table |TT1 we have reported core radii for each materials and configurations, deter- 
mined from anisotropic elasticity theory according to the definition given in [1]. It is usually 
assumed that an adequate core radius value for covalent materials is one quarter of the Burg- 
ers vector. Here, we got a good agreement in the sole case of silicon carbide. Otherwise, 
the calculated core radius is usually 10-15% (Ge and GaAs) to 30% (Si and diamond) larger, 
showing the limitations of such an hypothesis. 

Elements in the same column of the periodic table are expected to share a number of 
properties. Here, it is then interesting to compare energy differences, calculated for group IV 
materials. Results from our first principles simulations are reported in the Figure [2l as well 
as previous Tersoff calculations [25]. Going from Ge to diamond, the C core stability clearly 
increases, until C becomes the most stable configuration (diamond). A possible explanation 
is the presence of hybridized sp 2 atoms in the core. Such an electronic structure is highly 
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unfavorable for Ge and 4s and 4p orbitals compared to sp 3 (C is unstable in that case), while 
it greatly lowers the core energy for diamond and 2s and 2p orbitals. Comparison between 
first principles and semi-empirical results indicates that the variation obtained with the Tersoff 
potential is qualitatively correct. The sp 3 — ► sp 2 change in the electronic structure is then 
adequately modeled by this potential. However, the energy variation is too steep, C instead 
of A being obtained as the most stable configuration for SiC. AEab, calculated from first 
principles, increases going from Ge to diamond, and the B configuration is always unstable. 
The B core is characterized by broken covalcnt bonds along the dislocation line, and the 
energy variation simply suggests that it is more expansive to break bonds for diamond than 
for Ge. Using the Tersoff potential, B is metastable for all materials, and the energy variation 
is reversed. This discrepancy confirms that semi-empirical potentials are generally unable to 
describe with a good accuracy covalent systems including dangling bonds, such as surfaces or 
vacancy-cluster defects. 

In conclusion, by performing first principles calculations, we have shown that characteris- 
tics of a non-dissociated screw dislocation are very close for most of the investigated materials, 
i.e. silicon, cubic silicon carbide, germanium, and gallium arsenide. The most stable configu- 
ration is obtained for a dislocation lying in shuffle set planes. Only diamond exhibits different 
properties, the lowest energy structure being obtained for a screw dislocation located in the 
glide set. The stability of this glide core structure, including sp 2 atoms, improves significantly 
from Ge to diamond, a characteristic fairly well reproduced with the Tersoff potential. 

* * * 

One of us, L.P., wants to thank A. T. Blumcnau for fruitful discussions at the bar during 
the EDS2004 conference. 
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